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Abstract. Using first-principles calculations based on a variational density functional 
perturbation tlieory, we investigate the lattice dynamics of solid solutions of barium 
and strontium titanates. Averaging the information available for the related pure com- 
pounds yields results equivalent to those obtained within the virtual crystal approxi- 
mation, providing frequencies which are a good approximation to those computed for 
a (111) ordered supercell. Using the same averaging technique we report the evolution 
of the ferroelectric and antiferrodistortive instabilities with composition. 



INTRODUCTION 

Ferroelectric solid solutions with perovskite structure can exhibit properties 
which are much more interesting for technological applications than those of the 
related pure compounds. This is the case for the new class of single crystal re- 
laxor ferroelectrics, exhibiting exceptionally high piezoelectric constants [1]. It is 
also true for Baj^Sri.j-TiOs (BST), which combines a high dielectric constant with 
other properties making it a good candidate for the construction of dynamic ran- 
dom access memories (DRAM) [2]. 

In contrast to barium titanate (BaTiOs) which is a prototype ferroelectric (FE), 
strontium titanate (SrTiOa) is an incipient ferroelectric which undergoes an anti- 
ferrodistortive (AFD) phase transition at 105K [3] . Due to the different behavior of 
the related pure compounds, the structural properties of BST evolve with composi- 
tion [4]. It exhibits a sequence of ferroelectric transitions similar to that of BaTiOs 
for Ba concentrations ranging from 20 to 100%. Close to pure SrTiOs the behavior 
is complicated by the competition between FE and AFD instabilities yielding the 
phase diagram reported in Ref. [4]. 

In the past few years, first-principles effective Hamiltonians sucessfully described 
the temperature-dependent behavior of various pure perovskite oxides [5]. The 
generalization of this theoretical approach to solid solutions is a new challenge 



which includes some additional complications to treat the alloy disorder accurately. 

A straightforward procedure is to try to extract the relevant information on the 
alloy from computation for various ordered supercells. This idea was applied to 
the study of Pb^^Gei-j^Te solid solutions [6] but has the inconvenience of being 
computationally intensive. An alternative and promising approach is based on the 
Virtual Crystal Approximation (VGA) which allows the retention of a crystal with 
the primitive periodicity but composed of virtual atomic potentials averaging those 
of the atoms in the parent compounds. A careful construction of the VGA for 
PbZra.Ti(-i_2,)03 yielded a stress induced phase transition [7] and a compositional 
phase boundary [8]. Going further, the construction of a model Hamiltonian based 
on VGA calculations for the same compound recently was used to describe most of 
its temperature-dependent properties [9]. In the latter case, the VGA calculation 
provided a reference structure to which corrections had to be added to treat the 
behavior of the real ions correctly. 

If the purpose of the VGA is to define an average crystal from which the real 
solid solution can be described by inchiding a restricted number of small corrections, 
one can ask if a similarly good estimate of the mixed compound could be obtained 
directly from the information available for the pure materials by averaging some 
key quantities. This idea is reinforced by the observation that the interatomic force 
constants seem very similar in the different ABO3 compounds [10]. 

In what follows we define an Average Grystal Approximation (AGA) by averaging 
the lattice constant (ace//); ^^e Born effective charges (Z*), the optical dielectric 
tensor (eoo), and the interatomic force constants in real space (IPG) of the pure 
compounds. We compare our results to those obtained (i) within the VGA and 
(ii) for a cubic ordered supercell with alternating planes of Ba and Sr along the 
(111) direction. We show that our AGA yields an excellent approximation to the 
dielectric and dynamical properties of the mixed compound. It allows the prediction 
of the lattice dynamics in the full range of composition without requiring any first- 
principles calculations other than those for the related pure compounds. 

METHODOLOGY 

General framework 

Our calculations have been performed within the density functional theory [11] 
and the local density approximation, using the ABINIT package [12]. The ionic 
potentials screened by the core electrons have been replaced by highly transferable 
extended norm conserving pseudopotentials as proposed by M. Teter [13]. The 
Ba 5s, Ba 5p, Sr 4s, Sr 4p, Ti 3s, Ti 3p, Ti 3d, Ti 4s, O 2s and O 2p electrons 
have been treated as valence states. The wavefunction has been expanded in plane 
waves up to a kinetic energy cutoff of 45 Ha. The electronic exchange-correlation 
energy has been computed using a polynomial parametrization [14] of Geperley 
Alder data [15]. Integrals over the simple cubic Brillouin zone have been replaced 



by sums on a 6 x 6 x 6 mesh of special fc-points. A grid providing an equivalent 
sampling has been used for the supercell calculations. 

The Born effective charges, the dielectric constants and the dynamical matrices 
have been computed within a variational formulation of the density functional 
perturbation theory [16]. The interatomic force constants in real space (IFC) and 
the phonon dispersion curves have been extracted following the scheme described 
in Ref. [17], which deals correctly with the long range character of the dipolar 
interaction. The IFC were obtained from the knowledge of the dynamical matrix 
on a 2 X 2 X 2 BCC mesh of g-points. 



The virtual crystal approximation is based on the construction of a fictitious 
"virtual atom'' potential by averaging the ionic potentials of the atoms alternating 
at the same site of the structure. As illustrated recently by Ramer and Rappe [7], 
there are different ways of combining non-local atomic pseudopotentials, each of 
them leading to different physical results. A careful choice seems necessary to 
predict accurately the alloy properties at the VGA level. However, this choice could 
be less stringent if we think about the virtual crystal as a reference structure from 
which the real crystal can be obtained by including a limited number of relevant 
corrections. 

For EST, we choose therefore to stay with the most straightforward combina- 
tion of the reference atomic pseudopotentials. The potential at the A-site can be 
formally defined for a concentration of x% of barium as: 



where V^l and Vf^ are the pseudopotentials of Ba and Sr atoms. 

Non-local pseudopotentials, such as those we are using, can be expressed as a 
sum of a local contribution and short-range non-local corrections: 



Virtual atom pseudopotential 



V^'aA[x]-xVE: + {l-x)Vs': 



(1) 



(r)5(r-rO + EE 



l,m q 



Vr\4>lm,q >< <Plm,q\Vl 
< <f>lm,q\Vr\(t>lm,q > 



(2) 



Our virtual atom pseudopotential is then defined as follows [18]: 
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For X ranging from to 1, our virtual atom evolves smoothly from Sr to Ba. 



TABLE 1. Evolution of the LDA optimized lattice constant (bohr) 
with composition. Experimental values on the third line are interpo- 
lated from the data in the pure compounds using Vegard's law. 





SrTiOs 


x= 25% x= 50% 


x= 75% 


BaTiOa 


LDA/ VGA 


7.273 


7.321 7.370 


7.412 


7.451 


LDA/Supercell 




7.366 






Experiment 


7.38 


7.43 7.47 


7.52 


7.56 



GROUND-STATE PROPERTIES 



First we focus on the lattice constant of the cubic perovskite structure and look 
at its evolution with composition. In Table 1, we report the LDA lattice constants 
for BaTiOg and SrTiOs, for Bai.^Sr^TiOs within the VGA (for x = 0.25, 0.50 and 
0.75), and for the (111) ordered supercell. We compare the computed values with 
experimental results, which at intermediate compositions are accurately given by 
Vegard's law [4]. While our results show the lattice constant underestimate typical 
of LDA, the Vegard's law behavior is well reproduced. 

We note that the theoretical LDA underestimate of the experimental lattice 
constant fortuituously corresponds to the effect of a similar external pressure of 8.7 
GPa in both BaTiOa and SrTiOs. In what follows, we will work at the experimental 
lattice constants which is therefore equivalent to imposing a fixed fictitious negative 
pressure of -8.7 GPa to the LDA results. For the mixed compounds we will work 
at the values obtained by Vegard's law as reported on the third line of Table 1. 



DIELECTRIC AND DYNAMICAL PROPERTIES 
Born effective charges and dielectric tensor 

In Table 2, we compare the Born effective charges and the macroscopic optical 
dielectric constants obtained in our various calculations. The values of Oy and 0_l 
refer respectively to the change of polarization for a displacement of the oxygen 
atom along the Ti-0 direction and perpendicular to it. 

We observe that the values of Z* and eoo in the supercell are very similar to both 
those averaged from the pure compounds (AGA) and those obtained in the VGA. In 
particular, the AGA and the VGA similarly reproduce the large anomalous values 
of Z^j^ and which play an important role in generating the FE instability. For 
Z^, both kinds of approximation yield the average of the values computed for Ba 
and Sr in the supercell. A similar observation was reported by Bellaiche et al. for 
PZT [9]. 



TABLE 2. Born effective charges and optical dielectric con- 
stant in the cubic perovskite structure. In the supercell, we can 
have either Ba or Sr at the A site: the first value refers to Ba, 
the second to Sr. 
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Phonon frequencies and interatomic force constants 

Finally we compare the phonon dispersion curves, yielding information at the 
harmonic level on the full lattice dynamics of the crystal. Phonon frequencies are 
obtained by diagonalizing the dynamical matrix. The latter is usually decomposed 
into a dipolar part computed from the values of Z* and eoo, and a short-range part 
computed using the IFC. Our AC A dynamical matrix is constructed by simply 
averaging all the previous quantities from their values in the pure compounds and 
using an average mass at the A-site, while the VGA requires a full set of additional 
first-principles calculations with the virtual atom potential. 

In Fig. 1, we show the phonon dispersion curves of a mixed crystal of BST with 
x=50%, obtained both within the ACA and the VGA. We compare the results to 
those obtained for the related pure compounds. The dispersion curves of BaTiOs 
were discussed in Ref. [20]. Those of SrTiOa have been recomputed and agree with 
the calculation reported previously by LaSota tt al. [19]. Compared to BaTiOa 
which exhibits only a FE instabihty with a fiat dispersion in the F-X-M plane, 
the FE instability in SrTiOs is more localized in g-space (around the F point) and 
coexists with an AFD instability extending along the Brillouin zone edges (R-M 
lines) . 

We observe that averaging the different quantities (Z*, eoo, IFG, atomic masses) 
does not yield an average of the phonon frequencies but rather is closer to the 
results for pure BaTiOa. This is also what is found in the VGA, the two phonon 
dispersions being very similar. The major disccpancy appears for the low frequency 
Ti dominated branch along the F-M-X lines where the value of w^, being near zero, 
is especially sensitive to small changes in the dynamical matrix. 

In Table 3, we compare the AGA and VGA phonon frequencies at F and R in 
the reference perovskite structure to the related frequencies at the F point in the 
(111) ordered supercell. The different symmetry of the supercell does not allow a 
one to one correspondence but the comparison is made between phonons having 
the highest eigenvector overlap. We note that the supercell eigenmode associated 
with the Ri5 modes of the perovskite structure, and involving displacement of the 
A atom, exhibits a small LO-TO splitting. 




FIGURE 1. Phonon dispersion curves along high symmetry hnes of the simple cubic Brillouin 
zone for BaTiOa, Bao.sSro.sTiOs and SrTiOs. In the middle panel, the black lines correspond to 
the ACA, and the grey line to the VGA. 



TABLE 3. Phonon frequencies (cm~^) at the T and R points 
of the simple cubic Brillouin zone. Comparison of the ACA 

and VGA results with the supercell eigenmode at the F point 
which has the highest overlap with the corresponding ACA 
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ri5(T0) 


453 
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R25 


182 


84i 


135 
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R15 


414 
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766 


763 



The overall agreement between ACA, VGA and supercell calculations is very 
good. This is particularly true for the and R25' modes related to the FE 
instability. There is a larger discrepancy for the R25 AFD mode. We note however 
that projecting the supercell dynamical matrix within the subspace of the VGA 
eigenvector yields a value of 125 cm^^ (instead of 86 cm"^) closer to VGA and AGA 
results. This shows that the discrepancy originates in a change of the eigenvector 
rather than in a modification of the interatomic force constants. 



Evolution of the phonon instabilities with composition 



As the ACA seems to give accurate results for the lattice dynamics of BST sohd 
solutions, it offers an easy way to make a guess for the evolution of the FE and 
AFD instabilities with composition. We observe in Figure 2 that the square of the 
frequency of the R25 AFD mode evolves linearly with composition and becomes 
unstable for Ba concentrations smaller than 20 %. This corresponds precisely to 
the composition for which a change of behavior is observed experimentally in the 
phase diagram of BST [4]. In contrast, the FE instability exhibits a non-linear 
behavior and remains unstable for the full range of composition. 
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FIGURE 2. Evolution of the R25 and Fis unstable mode with composition, using the ACA. 

CONCLUSIONS 

In this paper, we have compared different approaches for the study of the lattice 
dynamics of BST sohd solutions and have introduced a new Average Crystal Ap- 
proximation. We have shown that the latter yields results equivalent to the VGA 
but with the advantage that it does not require any additional first-principles cal- 
culations if the information for the related pure compounds is available. Looking 
at the actual values of the IFC in the supercell could provide indications of the 
corrections to be included in the harmonic part of an effective Hamiltonian based 
on the ACA for the correct description of BST sohd solutions. 
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